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Abstract. Surge channels on wave-swept rocky shores 
are characterized by the violent hydrodynamic mixing 
that accompanies broken waves. It has been suggested 
that this mixing rapidly dilutes gametes shed into the surf 
zone, thereby severely reducing the fraction of eggs that 
can be fertilized externally. Although surge channels are 
well mixed within themselves, field experiments show that 
the exchange of water between these small embayments 
and the adjacent mainstream is surprisingly slow. Thus, 
surge channels may act as “containment vessels," limiting 
the rate at which gametes are diluted, and thereby en¬ 
hancing the efficacy of external fertilization. Indeed, a 
mathematical model of fertilization in surge channels 
suggests that given a sufficient population of adult males 
within a surge channel, 80-100% of eggs may be fertilized. 
This result must be tempered, however, by the possibility 
that the small-scale shears induced by turbulence interfere 
with fertilization. 

Introduction 

Recent field experiments with marine organisms have 
demonstrated that the fraction of eggs externally fertilized 
during spawning can be adversely affected by water mo¬ 
tion. For example, when water velocity is 0.1-0.2 m s 1 
and an individual male is separated from a female by 
more than 5 m, studies of sea urchins (Pennington, 1085; 
Levitan, 1991 ) and a hydroid (Yund, 1990) show that 
only 5—15 r of eggs arc fertilized. This low rate of fertil¬ 
ization is attributed to the dilution of gametes by turbulent 
mixing in flow . The percentage of eggs fertilized can be 
somewhat higher (approximately 30-50%) if multiple 
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males are present, but is lower at higher velocities (Pen¬ 
nington, 1985; Levitan et al, 1992). Even in fish that 
actively pair to spawn, the fraction of eggs fertilized can 
be low (40-50%) when the flow exceeds 0.1 m s^ 1 and 
wave-induced turbulence disperses the gametes (Petersen, 
1991; Petersen, et al. 1992). 

These studies (and the low rates of fertilization they 
report) raise important questions about the role of water 
motion in external fertilization and the consequences it 
may have in the evolution of reproductive strategies. 
However, flow during these experiments w'as relatively 
benign. How effective can external fertilization be under 
more extreme hydrodynamic conditions? 

As ocean waves break upon rocky shores they gen¬ 
erate water velocities as high as 15 ms ’( Denny. 1988; 
Vogel, 1981 ). twenty to 100-fold faster than those in 
the studies cited above. These high velocities are ac¬ 
companied by much more violent mixing of water. The 
“white water" characteristic of the surf zone is v isual 
ev idence of the extreme turbulence in breaking waves, 
and Denny and Shibata ( 1989) proposed that turbu¬ 
lence in the benthic boundary layer of the surf zone 
reduces the effectiveness of external fertilization well 
below even the low values measured subtidally. They 
estimate that only about 0.01-0.1% of eggs released into 
the surf zone are fertilized by a single adult male under 
typical conditions, and a maximum of only about 3% 
arc fertilized when multiple males spawn synchro¬ 
nously. These predicted rates of fertilization arc one to 
two orders of magnitude lower than those measured in 
subtidal Hows, and suggest even more dramatic con¬ 
sequences for intertidal plants and animals than those 
proposed for subtidal organisms. 

Denny and Shibata (1989) note, however, that the 
fraction of eggs fertilized in the surf zone could be sub- 
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stantiallv higher than they predict if the volume into which 
gametes can be diluted is limited. Experiments described 
here show that blind-ended surge channels (the small- 
scale embayments typical of rocky shores) can act to limit 
the dilution volume for gametes of benthic organisms. 
Shoreline topography may thereby enhance the effective¬ 
ness of external fertilization in the surf zone, and the frac¬ 
tion of eggs fertilized in surge channels may actually be 
similar to that found in subtidal habitats. 

Materials and Methods 

Field experiments were conducted on the shoreline ad¬ 
jacent to Hopkins Marine Station, Pacific Grove, Cali¬ 
fornia. Four blind-ended surge channels (essentially, 
small-scale bays) were chosen to span the size typical of 
this coastline. Each was surveyed with a surveyor's level, 
tape measure, and stadial rod to provide the location (in 
cylindrical coordinates) of 100-150 points, and these 
points were analyzed (Surfer, Golden Software, Inc.) to 
produce a topographic map of each surge channel (Fig. 
1). From these maps, the volume of each channel was 
estimated as a function of still water level (Fig. 2). Volume 
(at approximately the mean still-water level present during 
this study) ranged from 3.5 to 45 m 3 . 

The residence time of water in these channels was mea¬ 
sured at high tide as follows. A control sample of water 
in the center of the channel was taken with a test tube 
attached to a wooden pole. For this and all subsequent 
samples, the tube was filled and emptied once before sam¬ 
pling to minimize the effects of any ccntaminants in the 
tube. Thirty grams of water-soluble fluorescein dye was 
dissolved in a liter of seawater, and 100-200 ml of the 
solution was introduced into the surge channel. The dye 
was allowed to mix throughout the channel until it ap¬ 
peared uniform, a process resulting from the entry of one 
to three breaking waves into the channel (10-30 s). 
Dipped samples were then taken at intervals of 15 or 30 
s (depending on the rate of exchange between surge chan¬ 
nel and mainstream) until the dye was no longer visible 
in the channel. Typically 15-30 samples were taken in 
each trial. Two trials were conducted for each channel on 
each sampling day. Trials were alternated between sites 
to allow dye to be thoroughly dispersed from each site 
prior to the second trial. Sample tubes were capped and 
kept in the dark until analyzed (usually the next day). 

Markers were glued to one wall of each surge channel 
at 25 cm vertical intervals, allowing us to estimate visually 
(to within approximately 10 cm) the maximal and min¬ 
imal height of the water's surface during each surge. Mea¬ 
surements were recorded in each channel, for each surge, 
during a 5 min period both before the first trial and after 
the last trial. The average of each 5-min series of readings 
provided an estimate of the mean, or still-water, level in 
the channel during the trial. The average of minimal 


heights provided an estimate of the mean low water level. 
The number of surges encountered, divided by the 5-min 
duration of the recording session, provided an estimate 
of the average period of the surge. Periods ranged from 
about 9-14 s. The height of each surge was estimated as 
the maximal height minus the subsequent minimal height. 
Heights recorded in each 5-min period were averaged to 
provide an estimate, H avg , of the sea state during each 
trial. Average surge heights varied from 0.14 m to 
0.74 m. 

The concentration of dye in each sample w as measured 
with an Aminco model SPF 500 spectrofluorometer and 
a standard series of dilutions from the stock dye solution. 
Concentrations were corrected for any background fluo¬ 
rescence present in the control sample, and were expressed 
as a fraction, C r , of the concentration in the initial ex¬ 
perimental sample. A model of the time series of relative 
concentrations from each trial was then determined by 
fitting the data to the equation: 

C r (t) = a exp(-kt) + c; (Eq. 1) 

a simplex algorithm (Caceci and Cacheris, 1984) was used 
to estimate the least-squares fit. The parameter c represents 
the quasi-equilibrium concentration of dye at the end of 
each trial when the concentration in the surge channel is 
equal to that in the adjacent mainstream flow. In all cases, 
c is small (typically <0.000 la). The exchange parameter 
k represents the fraction of dye escaping from the surge 
channel in each interval of time t, and is thus a measure 
of the rate of exchange of water between surge channel 
and mainstream. The expected residence time, T r , of a 
dye particle in the surge channel is 

T r = 1/k. (Eq.2) 

A total of 73 trials were conducted between October, 
1988 and June, 1990. In three trials (4.1%), the measured 
concentrations in a surge channel appeared to increase 
substantially for one or more samples in the middle of 
the trial. These increases occurred too late in each trial 
to be attributed to incomplete initial mixing, and were 
due either to contamination in the sample tubes or to the 
accidental introduction of concentrated dye into the surge 
channel from an undetected “refuge" (a high tide pool, 
for example). Data from these anomalous trials were not 
included in the analysis. 

A mixing model for blind-ended surge channels 

As a means for examining our data, we formulated 
a simple model of the exchange between a surge channel 
and the mainstream. For the purposes of this model, 
the topography of a surge channel is represented by its 
volume versus still-water-level curve as shown in Figure 
2. At the mean low water level determined for a trial, 
the channel holds a volume V 0 . When a surge enters 
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Figure I. Topographical maps of the four embayments used in this study. Vertical gradations are 0.1 
m; hachures denote local depressions. 


the channel the water level rises, increasing the chan¬ 
nel's volume by an amount AV. If the initial volume 
(V 0 ) and added volume (AV) are thoroughly mixed, 
the probability that a particle is carried out of the chan¬ 
nel as the surge subsides and the water returns to its 
original level is 


Probability per surge = 


AV 

(V 0 + AV) * 


(Eq.3) 


But the mixing of initial and added volumes may not be 
complete; in such a case the probability of escape from 
the channel is only a fraction, m, (0 < m < 1) of this 
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shore that lie shoreward of this average were defined to 
be surge channels; segments seaward of the average shore¬ 
line were promontories. At each scale, the fraction of the 
average shoreline that comprised the mouths of surge 
channels. A, was then taken as a measure of relative surge 
channel abundance. This procedure was repeated for sev¬ 
eral starting points and the results averaged. The abun¬ 
dance of surge channels was measured at scales ranging 
from 2.2 to 22.4 m. 

The shore at Hopkins Marine Station is formed from 
granite, and is not atypical of rocky shores on the west 
coast of North America. 


Still Water Level (m) 

Figure 2. The relation between still water level and surge channel 
volume for the four embaymenls used in this study. These curves were 
calculated using the topographical maps of Figure 1. 


maximal probability. If the average period between surges 
is T seconds, the probability that a dye particle is carried 
out of the channel in one second is thus 


k = m 


AV 


(V 0 + AV)T 


(Eq.4) 


where k is again the exchange parameter of Eq. 1. 

V 0 . AV, T, and k can be determined experimentally, 
allowing us to calculate the mixing parameter m: 


m — 


k(V 0 + AV)T 
AV 


(Eq.5) 


Abundance of surge channels 

The abundance of surge channels on a shore is strongly 
dependent on the local topography, and therefore may 
vary greatly from site to site. Due to the fractal nature of 
shorelines (Mandelbrot, 1982: Denny, in press), it is dif¬ 
ficult to quantify precisely the abundance of surge chan¬ 
nels for even a single site, because the exact location of 
the shoreline varies with the scale at which the shore is 
examined. An estimate of the abundance of surge channels 
on the shores at Elopkins Marine Station was made, how¬ 
ever, using the following technique (Fig. 3). 

An aerial photograph of the shore was traced onto paper 
to provide a two-dimensional representation of the shore¬ 
line with accurate detail present at a scale smaller than 
that used in the analysis. A point was chosen at random 
on this trace, and a caliper placed with one tip at this 
point. The caliper was then “walked" down the trace, and 
the distance between the tips of the caliper set the scale 
at which the shoreline was measured. The scries of con¬ 
tiguous lines joining points where the caliper intersected 
the trace defined one realization of the “average" shoreline 
at that particular scale (Fig. 3). Segments of the actual 


Theory of external fertilization 

If the exchange between water in surge channels and 
that in the adjacent mainstream is slow, then surge chan¬ 
nels may function as “containment vessels" for gametes, 
thereby allowing for the effective fertilization of eggs. This 
possibility can be examined mathematically. 

In this analysis, we follow the model of Denny and 
Shibata (1989), based primarily on data from sea urchins. 
Males are assumed to release sperm into a surge channel 
at a constant rate Q s per second resulting in a concentra¬ 
tion of sperm C s (with units of m 3 ) that may vary through 
time. Denny and Shibata ( 1989) also allowed C s to vary 
in space, but for present purposes we assume that sperm 
are thoroughly mixed, ensuring that C s is the same every¬ 
where within the surge channel. Sperm (like dye particles) 
have a probability k in each second of passively escaping 
from the surge channel into the adjacent mainstream, and 
we assume that sperm that escape from the surge channel 
do not re-enter. 

Sperm can swim and, under certain conditions, can 
move preferentially toward eggs ( e.g Miller, 1985). The 
speed of sperm locomotion is very slow, however, a few 
tens of /um s ' 1 . In contrast, the motion of neutrally buoy- 



Figure 3. An example showing the traced shoreline (thick line), 
“average" shoreline at a particular scale (thin line segments), and the 
surge channel and promontory portions for a single segment of the average 
shore. This shoreline was traced from an aerial photograph of the shore 
at Hopkins Marine Station. 
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ant particles in turr \v can be characterized by a 

velocity u* (the aocity), a measure of the in¬ 
tensity of tu now (Schlichting, 1979; Denny, 

1988; Den /ata, 1989). For the intensely tur¬ 
bulent 11 i the surf zone of wave-swept shores, 

u* is ap; nel> equal to 0.1 u avg , where u avg is the 

ensemble rage velocity of the water (see Denny, 1988; 
Denny and Shibata, 1989). Water velocities in the surf 
zone vary from approximately 1-15 m s 1 (Denny, 1988), 
suggesting that u* is 0.1-1.5 m s' 1 . Thus, the velocity 
imposed by turbulence is four to five orders of magnitude 
faster than a sperm's swimming speed, and we assume 
that the swimming speed of sperm can be safely neglected 
in the present context. 

Given these assumptions, we see that the rate of change 
of sperm concentration is equal to the difference between 
the rate at which concentration increases due to the release 
of sperm and the rate at which concentration decreases 
due to the passive escape of sperm to the mainstream: 


dC 5 (t) 

dt 


y-kC s (t), 


(Eq.6) 


where V is the still-water volume of the surge channel. 

If we assume that C s (0) = 0, we calculate that 

QlD^^l 1 - exp(-kt)]. (Eq.7) 


At infinite time C s reaches a steadv-state concentration. 

C 5 . x : 

Q x = 777 • (Eq.8) 


We assume that females release eggs at the constant 
rate Q e per second, and that eggs are mixed and exchanged 
as for sperm. Thus, 


dC e (t) 

dt 


y~ kC e (t) 


(Eq.9) 


C e (t) = ^ [! - exp(-kt)], (Eq. 10) 


C e 


Qe 
kV ' 


(Eq. 11) 


where C e is the instantaneous concentration of eggs (m 3 ) 
and C c x is the steady-state concentration. 

We assume that the rate at which sperm fertilize eggs 
is governed by the co-occurring concentrations of both 
sperm and eggs and by a “reaction" parameter, 0. 

Rate of fertilization = 0C StV (t)C e . v (t). (Eq. 12) 


assumed to be the friction velocity u*) and the effective 
"collision cross section" of an egg (Denny and Shibata, 
1989), and has units m 3 s -1 . 

If an egg could be fertilized by the first sperm that con¬ 
tacted any point on its surface, the "collision cross section" 
would be the same as the projected area of the egg. For 
sea urchin eggs with a diameter of 10 4 m this is 7.85 
X 10“ 9 nr. Laboratory measurements suggest, however, 
that sea urchin eggs behave as if only a small fraction of 
their surface area is fertilizable. The reason for this reduced 
cross section is unclear. Vogel et al. (1982) suggest that 
apparent fertilizable area is only about 1% of the overall 
area for Pamcentwtus livid us. More recent experiments 
with Stwngylocentrotus fransiscanus (Levitan et al., 1991) 
suggest an apparent fraction of 3%. Using this latter value, 
we find that 

<t> = u*(2.4 X 10 10 m 2 ), (Eq. 13) 

and because u* 0. Iu avg , 

0 « u avg (2.4 X 10 11 nr). (Eq. 14) 

Given that average water velocities in the surf zone vary 
from approximately I to 15 m s 1 ,0 is likely to vary from 
2.4 X 10 11 to 3.5 X 10 10 m 3 s' 1 . 

Denny and Shibata (1989) calculated that an individual 
urchin extrudes sperm at approximately 10 s 1 or eggs 
at 10 4 s 1 . and these values are used here. Due to the vast 
overabundance of sperm it seems safe to assume that the 
concentration of sperm is minimally affected by the at¬ 
tachment of sperm to eggs; thus 

C SV = C S . (Eq. 15) 

Given these assumptions, we can describe the flux of virgin 
eggs through a surge channel. Virgin eggs are added to 
the local population as they are released by females and 
are removed either by escaping to the mainstream or by 
fertilization. Thus: 


dC r (t) = % - kCe. v (t) - 0C s (t)Ce. v (t). (Eq. 16) 
dt V 


At this point the mathematics can be simplified if we 
convert each of the variables of interest to a dimensionless 
form. In the present context, time is important primarily 
as it relates to the residence time of gametes in a surge 
channel. Thus we may specify a dimensionless time r, 

r = t/T r = kt. (Eq. 17) 

Note that 


dg = dg dt = I dg 
dr dt dr k dt 


(Eq. 18) 


Here, C s<v and C e<N are the concentrations of unattached 
sperm and virgin eggs, respectively. 

The reaction parameter is the product of the effective 
speed at which sperm are mixed through the water (here 


for any function g. 

The instantaneous concentration of eggs can be ex¬ 
pressed in dimensionless form by dividing by the steady- 
state concentration of eggs. Thus, 
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C e (r)_ C c (r) 
Ce.x (Qe/kV) ' 


(Eq. 19) 


Similarly, the instantaneous concentration of virgin eggs 
is: 


C, v (r) C e v ( r ) 
Ce.x (Qe/kV) ' 


(Eq.20) 


The instantaneous concentration of sperm is made di¬ 
mensionless by expressing it as a fraction of the steady- 
state concentration of sperm: 


_ C s (t) _ C s (t) 

= Cs.* (Qs/kV) ' 


(Eq. 21) 


Finally, we devise a dimensionless fertilization parameter 
0 that incorporates several of the variables that govern 
the rate of fertilization: 


" k sx k : V V 


(Eq.22) 


Inserting these dimensionless variables into our previous 
results (Eqs. 7, 8, 10, 1 1. and 16), we find that: 


s(r) = e( t ) = 1 - exp(-r), (Eq.23) 
s x = e* = 1, (Eq.24) 

= 1 - v(t)-9s(t)v(t). (Eq.25) 

dr 


Thus, at steady state (that is, when dv(r)/dr = 0), 

v* = l/d +©)• (Eq.26) 

In other words, the larger 0 is, the lower the steady-state 
fraction of virgin eggs, v x . A large 0 can result from a 
large rate of sperm release (for example, due to the pres¬ 
ence of multiple males), a small surge channel volume, 
an intense level of turbulence, or a long residence time 
(Eq. 22). Thus, any combination of these factors produces 
a lowered fraction of virgin eggs. Note, however, that 0 
scales with the square of residence time, so v x is especially 
sensitive to T r (or, equivalently, to k). 

Now, the concentration of fertilized eggs is equal to the 
difference between the total concentration of eggs and the 
concentration of virgin eggs. By analogy, the fraction of 
all eggs that are fertilized, f(r), is equal to 1 — v(r). From 
Equation 26 we see that at steady-state 

f, -0/(1 +0). (Eq.27) 

Xon-steady-state conditions 

This result applies to surge channels once they have 
reached a steady-state concentration of gametes. Is it realistic 
to apply these steady-state results to the real world? To an¬ 
swer this question, we solved Equation 25 for v as a function 
of r using an implicit trapazoidal finite difference technique: 


v([i + 1]3r) 

{1 - (5r/2)[ l + 0 ~ 0exp - (i hr + ^r)]} 

” V( ' 5t) {1 + (<5r/2)[l + 0 - 0exp - (i5r + Vibr)]} 

<5r 

+ {1 + (5r/2)[l + 0 - 0exp - (i Hr + */ 2 5r)]} ‘ 

(Eq.28) 

where hr is the time step and i is a counter (i = 0, 1, 2, 
3,. . .). At r = 0, all eggs are assumed to be virgin so that 
v(0) = 1. 

We again set f = 1 — v and follow fas a function of 
dimensionless time. 


Results 

Mixing in surge channels 

The decrease in dye concentration within a surge channel 
is modeled accurately by the exponential decrease of Equa¬ 
tion 1 (Fig. 4). The average coefficient of determination 
(r 2 ) for the 70 trials is 0.955 (with a standard deviation of 
0.054). indicating that, to a first approximation, each dye 
particle can be regarded as having a fixed probability of 
escaping from the surge channel in a given interval. 

Exchange between surge channel and mainstream water 
is surprisingly slow. The exchange parameter k varies from 
approximately 0.001 to 0.04, implying average residence 
times, T r , of 25 to 1000 s. 

The mixing parameter, m. varies from about 0.05 to 
0.9. When data from all four surge channels are combined, 
m is found to be an increasing function of average surge 
height (Fig. 5), but variation in H avg explains only 17% 
of the variation in m. The mixing parameter is not sig¬ 
nificantly correlated with the average fractional change in 
volume during a dowmsurge, AV/(V 0 + AV) (r = 0.0619, 
68 df, P> 0.5). 

When data from all four surge channels are combined, 
the exchange parameter k is found to be an increasing 
function of average surge height (Fig. 6), and the rela¬ 
tionship explains 51% of the variation in k. The exchange 
parameter is also positively correlated with the fractional 
change in volume in the surge channel, which is itself a 
direct function of H avg , but the correlation (k = 0.179 
H avg + 0.0023, r = 0.431, 68 df) explains much less of 
the variation than does surge height alone (19 vx. 51%). 

Abundance of surge channels 

Surge channels form about 46% of the shoreline at 
Hopkins Marine Station (Table 1). The fractional abun¬ 
dance of surge channels is not correlated with the scale 
at which the shoreline is measured (r = 0.134, 3 df, P 
> 0.5). at least for the range of scales examined here. 

Predicted fraction of eggs fertilized 

At a realistic value of 0 (that is, a realistic combination 
of k, V, </>, and Q s ), a substantial fraction of eggs is expected 
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Figure 4. Representative examples of the temporal decay in dye concentration in the experimental 
embaymenls. The solid curves are least-squares fils to the data using the model of Equation 1. 


to be fertilized (Fig. 7). For example, if 0 is 10 10 s~ l , 
and the residence time is 50 s (k = 0.02 s _1 ) for a surge 
channel with a volume of 10 m 3 (typical values for the shore 
at Hopkins Marine Station), about 20% of eggs can be fer¬ 
tilized by a single male. This is two to three orders of mag¬ 
nitude larger than the fraction predicted by Denny and Shi- 
bata (1989) for flow outside of surge channels. 

Multiple males 

Increasing the number of males present in a surge channel 
increases Q s , with a concomitant increase in 0, and thereby 
an increase in the fraction of eggs fertilized (Fig. 8). Consider, 
for instance, a surge channel with the same residence time 
as that cited above (k = 0.02 s ‘), but with five times the 
volume (50 m 3 ). If this channel has an average still water 
depth of m, it has approximately 70 nr of substratum 
area. Even ven the fivefold increase in volume, if only 3 
males are p nt per square meter of substratum (210 males 
total, a realistic \alue for wave-swept shores), 90% of eggs 
are expected to be fertilized. 

Smaller channels w ill have a larger ratio of substratum 
area to channel volume. Thus, given the same spatial 


density of male urchins, smaller channels will have a 
higher fraction of eggs fertilized. 

Time dependence of fertilization 

The time-dependent fraction of eggs fertilized is shown 
in Figure 9a. For realistic values of0, values off within 10% 
of steady-state values are obtained within a period of 5 T r . 
For a typical residence time of 50-120 s, this implies that 
the steady-state estimate of fertilization fraction is reached 
within 4-10 minutes. Given that an individual urchin in 
the laboratory releases gametes for approximately 1 h, equi¬ 
librium conditions are likely to apply to most gametes re¬ 
leased. 

Note that Figure 9a show's the fraction of eggs fertilized 
relative to the steady-state fraction. The absolute fraction 
fertilized is shown as a function of dimensionless time in 
Figure 9b. The fraction fertilized decreases with decreasing 
0 as expected. 

Discussion 

Mixing in small surge channels 

The rate of exchange between small well-mixed surge 
channels and mainstream waters is both a simple and a 
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Figure 5. The relation between average surge height (H avg ) and the 
dimensionless mixing parameter, m. The linear regression (solid line) 
explains approximately 17% of the variation in nr. m = 0.424 X H avg 
+ 0.094 (r 2 = 0.173), where H avg is measured in m. 

complex process. Our data indicate that the probability per 
time that a small, inanimate, neutrally buoyant particle will 
escape from a surge channel is constant for at least a few 
tens of minutes, and that this probability can be predicted 
with reasonable accuracy from the height of the surge present 
in the surge channel. In this respect the results of the process 
are simple to predict as shown in Figure 6. 

It is more difficult to account for these predictions on 
a mechanistic basis, primarily because of the variability 
of the mixing parameter, m. At present, m cannot be pre¬ 
dicted accurately, from either average surge height (Fig. 
4) or fractional change in volume. Probably, m is affected 
by the local topography of each surge channel and by 
chance variation in the shape of each breaking wave, ren¬ 
dering it of little value as a general predictive tool. 

Note that the relationship measured here between surge 
height and exchange parameter (Fig. 6) can be applied 



Figure 6. The relation between average surge height ( H avg ) and the 
exchange parameter k. The linear regression (solid line) explains ap¬ 
proximately 51% of the variation in k: k = 0.036 X H avg + 0.002 (r 2 
= 0.514), where H avg is measured in m, and k has the units s' L 


Table 1 


Fraction of the average shoreline that comprises the mouths of surge 
channels at Hopkins Marine Station 


Scale (m) 

Fraction, A 

2.2 

0.446 

4.5 

0.440 

9.0 

0.508 

22.4 

0.455 


only to a relatively narrow range of channel sizes. There 
is a lower limit to the size at which a local indentation in 
the shoreline can be considered a surge channel for our 
purposes. If the still-water depth in the indentation is less 
than the amplitude of the surge (=H avg /2), the indenta¬ 
tion is completely emptied by an average downsurge, and 
the probability of a particle re-entering the channel on 
the subsequent upsurge is set more by mainstream mixing 
and advection than by mixing within the indentation itself. 
In this case, the indentation is not a surge channel in the 
sense that we have used the term here. 

The upper limit to the size of this relationship is set by 
the size at which a surge channel is no longer well mixed 
within itself. Only if the channel is well mixed can Equa¬ 
tion 1 accurately model residence time. We have noticed 
distinct advective patterns in bays in excess of the size 
studied here—particles that enter a larger bay at one por¬ 
tion of its mouth predictably exit at another specific re¬ 
gion. In this case, the probability of exchange for a particle 
depends in part on its history, and the simple model of 
Equation 1 cannot be expected to hold. Thus, the pre¬ 
dictions implicit in Figure 6 should not be applied to surge 
channels larger than approximately 100 m 3 . 



k 

Figure 7. As the probability of gametes escaping from a surge channel 
increases, ihe equilibrium fraction of eggs fertilized decreases. 1 he larger 
the volume of the surge channel (V), the lower the fraction of eggs 
fertilized. 
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Figure 8. The fraction of eggs fertilized increases with the number 
of males present in a surge channel. 


The fraction of eggs fertilized 

The predictions made here suggest that external fertil¬ 
ization in surge channels may be quite effective, with fer¬ 
tilization fractions comparable to or exceeding those 
measured for benthic invertebrates in relatively benign 
subtidal flows (Pennington, 1985; Yund, 1990; Levitan, 
1991; Levitan et ai , 1992). If these predictions are valid, 
fertilization may be much less of a limiting process than 
suggested by the model of Denny and Shibata (1989), at 
least within surge channels. 

Furthermore, about 46% of the shoreline at I lopkins Ma¬ 
rine Station is formed of surge channels (Table II). If we 
assume that the fraction of eggs fertilized within the surf 
zone, but outside of surge channels, is equal to the maximal 
fraction predicted by Denny and Shibata ( 1989) (3%), and 
that individuals are randomly distributed along the shore, 
then the overall fraction of eggs fertilized is approximately 

Overall fraction = (f — 0.03)A + 0.03, (Eq.29) 


and sperm are mixed, dilution of gametes is so rapid that 
those few' eggs that are fertilized are likely to be fertilized 
by sperm from the nearest males. In this fashion, rapid 
dilution insures the localization of fertilization. 

Asa corollary 7 to the effectiveness of fertilization in surge 
channels, we note the implication that the long residence 
time of water in channels can again lead to localization 
of fertilization. Once gametes escape from a channel they 
are subject to the infinite dilution described by Denny 
and Shibata (1989), and subsequent fertilization becomes 
highly unlikely. Eggs released in a surge channel are 
therefore very likely to be fertilized by males in the same 
channel. Thus both infinite dilution and the limitation of 
dilution result in localized fertilization, leading us to sup¬ 
pose that, as a general rule, external fertilization in the 
surf zone can be effective only over short distances. 

The effect of localized fertilization due to surge channels 
on the population genetics of a species is unclear. Although 
the residence time of small surge channels is surprisingly 
long, it is still short compared to the larval lifetime of 
most species, which vary from a few hours to several 



where f is again the steady-state fraction of eggs fertilized 
in surge channels, and A is the fraction of shoreline that 
comprises the mouths of channels. For f = 0.9 and A 
= 0.46, the overall fraction of eggs fertilized is thus 43%, 
suggesting that external fertilization in the surf zone, taken 
as a whole, may be considerably higher than that suggested 
by Denny and Shibata (1989). 

The role of shoreline topography in controlling the overall 
rate of external fertilization is clear in Equation 29. If the 
shore is constructed of rocks that are not susceptible to the 
formation of surge channels, A will be small, and the overall 
fraction of eggs fertilized will be similarly small. Conversely, 
on a shore where surge channels are abundant, external fer¬ 
tilization ma be quite effective. 

Local i zed fertil i zat ion 

Denny and Shibata (1989) concluded that, as long as 
there is no practical limit to the volume into which eggs 



Figure 9. (A) Steady-state conditions are approached within a period 

equal to 5 limes the retention lime(r = 5). The fraction of eggs fertilized 
is expressed as a fraction of the steady-state value. (B) The actual fraction 
of eggs fertilized is plotted without being normalized to the fraction at 
steady state. 
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months. As a result, although fertilization may be local¬ 
ized, the resulting larvae are free to participate in the 
larger-scale dynamics of the population. We recognize that 
the long residence time of surge channels may have im¬ 
portant implications for the spatial pattern of settlement 
of larvae, as well as for the probability that newly fertilized 
embry os will be consumed by benthic suspension feeders, 
but these implications will not be developed here. 

Instantaneous mixing, shear stress, and a note of 
caution 

We foresee two factors that could affect the results pre¬ 
sented here. 

First, the model of external fertilization assumes that 
gametes are instantaneously mixed throughout a surge 
channel. As noted, however, thorough mixing within a 
channel requires the entry of one to three breaking waves. 
In the brief period prior to the arrival of the first wave, 
local high concentrations of gametes may be present, and 
if eggs encounter an area of high sperm concentration, 
the probability of fertilization may be very high. Note, 
however, that any momentary' reduction in mixing that 
allows a high concentration of sperm to be maintained 
also reduces the probability that an egg will be delivered, 
by chance, to that area of high concentration. Quantitative 
analysis of the dynamics of fertilization over times shorter 
than the period of the surge will require a more complete 
understanding of the structure of turbulence within surge 
channels than is currently available. 

Second, the calculations made here are based on the 
assumption that sperm and eggs interact in a fashion akin 
to molecules in a gas. If a sperm happens to be in the 
vicinity of an egg and aimed in the right direction, it im¬ 
pacts on the egg, sticks, and (perhaps) fertilization follows. 
However, this simple view ignores the small-scale hydro¬ 
dynamics relevant to sperm and eggs, and as a result may 
overestimate the probability of fertilization. 

To see why, we explore flow in the surf zone at the 
spatial scale of a gamete. This exploration is carried out 
in two steps. First, we examine the rate at which turbulent 
kinetic energy is dissipated in the surf zone. Next, we use 
these results to estimate the average velocity gradient to 
which gametes are subjected. The biological consequences 
of this gradient can then be discussed. 

Ocean waves carry' energy with them as they approach 
the shore. The rate, P, at which wave energy is transported 
across the mouth of a surge channel M meters wide is 

(1 /8)pgH 2 M(2gH) 1/2 

« 0.177pg 3/2 H 5/2 M, (Eq.30) 

where p is the density of seawater (approximately 1025 
kg m 3 ), g is the acceleration due to gravity (9.81 m s _l ), 
and H is wave height (Denny, 1988). Thus, the power 
per meter of surge channel mouth, P', is 


P' « 0.177pg 3/2 H 5/2 . (Eq.31) 

In making this calculation we have assumed that wave 
height at breaking is equal to water depth (Denny, 1988), 
and that waves at breaking move at the speed predicted 
by solitary' wave theory (in this case, (2gH) i; " [Denny, 
1988]). 

Some of this incident wave energy may be carried back 
out of the surge channel if the wave reflects from the shore. 
If the height of the reflected wave is a fraction R of the 
incoming wave height (0 < R < 1), the net power deliv¬ 
ered to the channel (per meter of channel mouth) is 

P' n » 0.177pg 3/2 H 5/2 (l -R 5/2 ). (Eq.32) 

On rocky shores, R (the refection coefficient) is typically 
small. For example, the R for a rubble breakwater (similar 
in topography to a rocky shore) is less than 0.29 when 
the slope of the breakwater is 1:5 (U.S. Army Corps of 
Engineers, 1984). If the slope is 1:10, the reflection coef¬ 
ficient is less than 0.13. If energy is transported inshore 
by waves but only a small fraction is carried back out by 
waves, most of the incident wave power must either be 
dissipated in the surf zone or transported out by other 
means (e.g., kinetic energy associated w’ith an undertow, 
acoustic energy). Effective transport by these other mech¬ 
anisms appears unlikely, and we assume here that most 
wave energy is dissipated in the surf zone. 

The magnitude of this rate of energy dissipation can be 
estimated from a simple example. The bottom of a surge 
channel has a slope a , as shown in Figure 10, and the wave 
height at the mouth of the channel is equal to the depth at 
the mouth of the channel. Given these conditions, waves 
break near the channel mouth, and the entire net energy' of 
the waves is dissipated within the channel. 

The still-water volume of the surge channel per w idth 
of mouth is 


As a breaking wave moves into the channel, its height 
decreases as wave energy is first converted to tubulent 
kinetic energy and then dissipated. Thornton and Guza 
(1983) have shown that the average (root mean square) 


Wave Height, H 



Figure 10. A schematic cross section through a surge channel illus¬ 
trating the relation between bottom slope (a) and water depth, and 
showing the typical shape of a broken wave. 
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height of waves afte- rreak on a given shore is a 
constant fraction a water depth. For a sandy 

beach, this fra- < 1 and it is likely to be similar 

on rock\ s' ... e of computation, however, we 

use a v . >> the broken wave propagates into 

the ch' mme that its shape approximates that 

of a uiri ore; that is, it has a steep leading face, and 

i level behind the face is nearly equal to that of 
the wave crest (Fig. 10). Given these assumptions, the 
broken wave increases the volume of water in the surge 
channel by maximally the channel’s still-water value. 
Thus, the maximal effective volume of water in which 
the wave energy is dissipated is 


V' 


5H 2 
8 tan a 


(Eq.34) 


The minimal rate at which wave energy is dissipated per 
unit of effective water volume (W') in the surge channel 
is thus, 

W' = P'J V' « 0.283pg 3/2 H ,/2 ( 1 - R 5/2 ) tan a. 

(Eq. 35) 


o = pft. (Eq. 37) 

If the cube has sides of length dx, the force required to 
shear the cube is thus 

F = dx 2 p/3. (Eq.38) 

The power, P, dissipated by this force is equal to the prod¬ 
uct of force and velocity. The velocity at the top of the 
cube differs from that at the bottom by dxft. Therefore, 

P = d xW 2 - (Eq. 39) 

The power dissipated per volume (dx 3 ) is thus 

W' = p/? 2 , (Eq.40) 

and the power dissipated per mass is 

( = pft 2 /p. (Eq.41) 

Equating Equations 35 and 40 and solving for /3, we find 
that the average velocity gradient in the fluid in a surge 
channel is: 

/3 avg « 0.532p 1/2 g 3/4 H l/4 ( 1 - R 5/2 ) l/2 (tan o:) l/2 P ]/2 < 

(Eq.42) 


Now, the turbulent kinetic energy due to wave breaking 
is ultimately dissipated as heat via the action of viscosity. 
This mechanism is elegantly described by Lazier and 
Mann (1989): at small spatial scales, turbulent eddies in 
the broken wave are damped by the water's viscosity, 
forming in their stead a linear velocity gradient (or shear) 
within the fluid that varies randomly in direction and 
strength. It is the interaction of this shear with viscosity 
that converts turbulent kinetic energy to heat. 

The scale at which the chaotic motion of turbulent ed¬ 
dies is effectively damped can be estimated from the Kol¬ 
mogorov length scale, L v : 

Lv = (p^) / 4 ' (Eq. 36) 

where p is the dynamic viscosity of seawater (1 to 2 
X 10 -3 Pa s, depending on temperature), p is again the 
water s density, and e is the rate at which turbulent kinetic 
energy is dissipated per mass of fluid (W kg -1 ). The largest 
shears are found in eddies with a diameter equal to approx¬ 
imately 40 Ly. Smaller eddies are much less energetic, and 
few eddies have a diameter less than 5-10 Lv (Lazier and 
Mann, 1989). 

The fact that wave energy is dissipated via small-scale 
velocity gradients provides a method for examining the 
characteristics of flow on the spatial scale of eggs and sperm. 
Consider a small cubic volume of water with its top face 
moving relative to its bottom face such that a linear velocity 
gradient ft (with units ofs 1 ) is established within the cube. 
As a result, the water in the cube is sheared. From Newton's 
law of viscosity (Vogel, 1981), we know that the shear stress 
(force per area) required to deform water at this rate is 


and (from Eq. 37) the average shear stress in the fluid is 

c7 avg « 0.532p ,/2 g 3/4 H 1/4 (l - R 5/2 ) l/2 (tan a) ,/2 n' /2 . 

(Eq. 43) 

Note that these are average values for shear and shear 
stress. At any point in the fluid, instantaneous values likely 
vary through time. For each interval where the shear stress 
is lower than average, however, there must be a compen¬ 
sating interval in which stress is larger than average. 

We are now in a position to evaluate directly the small- 
scale flow regime in surge channels and to explore its bio¬ 
logical consequences. 

For a slope of 1:5, a breaking wave height of 1 m, a re¬ 
flection coefficient of 0.29, and a dynamic viscosity of 1 X 
10" 3 Pa s, the average rate of energy dissipation within a 
surge channel is 1.66 W kg 1 (Eqs. 40 and 41). From this 
value we estimate that the eddies associated with maximal 
shear stresses are approximately 1094 pm in diameter, 11- 
14 times the diameter of sea urchin eggs (80-100 pm). The 
smallest eddies likely to be found in surf-zone flow are about 
137-273 pm in diameter. Even these small eddies are larger 
than typical gametes, suggesting that gametes in the surf 
zone experience a flow regime characterized primarily by 
linear shears rather than by turbulent eddies. 

Note that the heat capacity of water is about 4200 J 
kg" 1 K _1 (Weast, 1977). At the dissipation rate in this 
example, water is heated by only 0.0004°C s' 1 . Thus, the 
energy dissipation within a surge channel is not sufficient 
to raise the temperature substantially, even for retention 
times of 100-1000 s. 

For the typical conditions cited above, the average shear 
stress in the water is about 1.30 Pa, and the average mi- 
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croscalc velocity gradient is 1304 s 1 . The presence of 
this gradient means that the velocity at one side of an egg 
is quite different from that at the opposite side, which 
causes the egg to rotate. The rate of rotation (in cycles 
per second) is (Happel and Brenner, 1983; Kessler, 1986) 

T-. (Eq.44) 

4 7T 

In other words, an egg in a typical surf zone rotates on 
average about 104 times per second! Sperm oriented parallel 
to the velocity gradient (/>., the velocity at the tip of their 
flagellum is maximally different from that at their head ) will 
also be rotated by the shear, but because sperm are elongated 
rather than spherical, rotation is likely to cease once the 
sperm lies perpendicular to the instantaneous gradient ( i.e., 
with its head and tail in flow' at the same speed). 

The practical results of the rapid rotation of eggs, the 
rapid re-orientation of sperm, and the large shear stresses 
in the fluid are difficult to estimate. Certainly the egg, as 
a solid object, will entrain fluid with it as it rotates, creating 
a complex, 3-dimensional, temporally variable boundary 
layer that must be crossed by sperm. Whether the rotation 
itself can affect contact between egg and sperm is less clear. 
Also unclear is whether surf-zone shear stresses are suf¬ 
ficient to dislodge sperm from the egg once they have 
attached, or to damage the sperm or egg directly, although 
Denny and Shibata (1989) note preliminary evidence that 
even low shear stresses can reduce the fraction of eggs 
fertilized in the laboratory'. 

While the precise effects of the small-scale flow regime 
are unclear, the motions of sperm and egg will surely be 
complex and highly dynamic in turbulent flow which may 
also impede effective contact and subsequent fertilization 
and development. As long as the effects of shear stress 
and egg rotation remain unknown, the calculations made 
here concerning the effectiveness of fertilization in surge 
channels must be taken with a grain of salt. If turbulence 
inhibits the contact and subsequent attachment of sperm 
to eggs, or if the rapid rotation of eggs is detrimental to 
their survival, the calculations made here may overesti¬ 
mate the effectiveness of fertilization. 
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Symbol 


A 

C e (t; 
C, , < ! 

C,., 
c r (t) 
c s (t) 

Q, 

Q* 

e(r) 

e* 

f(r) 

f* 

F 

g 

H 

k 

U 

m 

P 

P' 

P', 

Qe 

Qs 

R 

S (r) 

S x 

t 

T 

T r 

U av g 

U* 

V (r) 
Voc 

V 
Vo 
V' 

v; wl 

A V 
W' 


a 

(3 

/^avg 

€ 

0 

<T 

^avg 

T 

4> 


s used in the text and the equation where each is defined or first used 


Definition 

Units 

Equation 

Relative abundance of channels 


29 

Concentration of eggs at time t 

m 3 

9 

Concentration of virgin eggs 

m“ 3 

16 

Steady-state egg concentration 

m 3 

11 

Relative dye concentration 


1 

Concentration of sperm at time t 

m 3 

6 

Concentration of unattached sperm 

m 3 

15 

Steady-state sperm concentration 

m 3 

8 

Dimensionless egg concentration 


19 

Steady-state e 


24 

Fraction of eggs fertilized 

Steady-state fraction of fertilized eggs 


27 

Force 

N 

38 

Acceleration due to gravity 

m s 2 

30 

Surge height 

m 

31 

Exchange parameter 

s _1 

1 

Kolmogorov length scale 

m 

35 

Mixing parameter 


4 

Energy flux into a surge channel 

W 

30 

Energy flux per meter of channel mouth 

W m 1 

31 

Net P' 

W m 1 

32 

Rate of egg release 

s 1 

9 

Rate of sperm release 

s 1 

6 

Reflection coefficient 


32 

Dimensionless sperm concentration 


21 

Equilibrium s 


24 

Time 

s 

1 

Period of the surge 

s 

4 

Expected residence time (1/k) 

s 

2 

Mean water velocity 

m s _l 

14 

Friction velocity 

m s 1 

13 

Dimensionless concentration of virgin eggs 


20 

Equilibrium v 


26 

Channel volume at still water level 

m 3 

3 

Surge channel volume at mean low water 

nr 1 

3 

Maximal volume per meter of channel mouth 

nr 

34 

Volume per meter channel mouth at SWL 

nr 

33 

Surge-induced increase in channel volume 

m 3 

3 

P'„/V' 

W m~ 3 

35 

Bottom slope 


33 

Velocity gradient 

s 1 

37 

Average velocity gradient 

s 1 

42 

Rate of energy dissipation 

Wkg 1 

36 

Dimensionless interaction parameter 


22 

Dynamic viscosity of water 

Pas 

36 

Water density 

kg m 3 

30 

Shear stress 

Pa 

37 

Average shear stress 

Pa 

43 

Dimensionless time 


17 

Reaction parameter 

m 3 s' 1 

12 
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